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PREFACE 
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height  given  the  straight-line  fetch  length,  the  vindapeed,  and  the  nearshore 
water  depth.  The  wave  height  prediction  curves  were  generated  by  numerically 
propagating  offshore  JONSWAP  spectra  shoreward  while  applying  shoaling  and 
wave  steepness  limitation  criteria  to  each  spectral  component.  The  report  pro¬ 
vides  an  alternate  approach  to  the  problem  of  shallow-water  wave  estimation. 

The  work  was  carried  out  under  the  shallow-water  wave  transformation  program 
of  the  U.S.  Army  Coastal  Engineering  Research  Center  (CERC) . 

The  report  was  written  by  William  G.  Grosskopf,  Hydraulic  Engineer,  and 
Dr.  C.  Linwood  Vincent,  Chief,  Coastal  Oceanography  Branch,  Research  Division. 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 

U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x  10"3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

lTo  obtain  Celsius  (C) 

temperature  readings 

from  Fahrenheit  (F)  readings, 

use  formula:  C  -  (5/9) 

(F  -32). 

To  obtain  Kelvin  (K) 

readings,  use  formula: 

K  -  (5/9)  (F  -32)  +  273.15. 

SYMBOLS  AND  DEFINITIONS 


Cfc  wave  steepness  limitation  factor 
d  water  depth 
E  energy  density 

EX  total  energy  in  the  wave  spectrum 

F  straight-line  fetch  length  (for  irregularly  shaped  water  bodies,  this 
should  be  based  on  an  average  over  a  24°  quadrant) 

f  frequency  of  spectral  component 

fm  frequency  of  spectral  peak 

g  acceleration  due  to  gravity 

Ho  deepwater  significant  wave  height 

Hs  significant  wave  height 

H  dimensionless  wave  height 

Ka  shoaling  coefficient 

L  wavelength 

R  land-wat~r  windspeed  correction  factor 

RT  air-sea  temperature  difference  windspeed  correction  factor 
Tp  peak  wave  period 

UA  windspeed  to  be  used  in  wave  height  estimation 
U&  overwater  windspeed  corrected  for  wind  instabilities 
Uw  overwater  windspeed 

Uz  windspeed  measured  Z  meters  above  land  or  water  surface 
U^o  10-meter  (33  foot)  windspeed 

x  dimensionless  fetch  length 

a  Phillips  equilibrium  constant 

Y  ratio  of  maximal  spectral  energy  to  the  maximum  of  the  corresponding 
Pierson-Moskowitz  spectrum 

it  3.14159 

left-side  width  of  the  spectral  peak 
Cfb  right-side  width  of  the  spectral  peak 
<j>  wave  steepness  limitation  factor 

<i>h  wave  steepness  limitation  factor 
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ENERGY  LOSSES  OF  WAVES  IN  SHALLOW  WATER 


„  .  by 

William  G.  Grosskopf  and  C.  Linuood  Vincent 
I .  INTRODUCTION 

The  energy  in  an  irregular  wave  train  changes  as  the  waves  propagate  from 
deep  water  toward  shore.  Estimates  of  the  total  change  in  wave  energy  have 
traditionally  been  made  by  multiplying  a  shoaling,  refraction  and  friction 
coefficient  by  an  offshore  significant  wave  height  to  yield  the  nearshore  wave 
height.  Recent  studies  of  wave  spectra  have  provided  a  more  detailed  view  of 
the  wave  field  and  Indicate  that  additional  processes  should  be  considered. 

This  report  presents  finite-depth  wave  height  estimation  curves,  given  an  ini¬ 
tial  JONSWAP  type  of  offshore  spectral  wave  condition  (Hasselmann,  et  al,,  1973) 
generated  over  a  short  fetch  and  incorporating  finite-depth  steepness  effects 
based  on  a  study  by  Kitaigorodskii,  Krasitskii,  and  Zaslavskil  (1975).  These 
curves  represent  energy  changes  due  to  shoaling  and  an  upper  limit  of  energy 
spectral  density  as  a  function  of  wave  frequency  and  water  depth. 

Research  at  the  Coastal  Engineering  Research  Center  (CERC)  and  elsewhere 
indicates  steepness  effects  that  lead  to  breaking  in  a  shoaling  wave  field 
lead  to  a  major  loss  of  energy  in  addition  to  that  lost  by  bottom  friction  and 
percolation.  These  effects  can  be  incorporated  into  wave  estimation  curves  in 
a  fashion  similar  to  shoaling  because  the  effects  can  be  made  a  function  of 
depth.  The  effects  of  refraction,  bottom  friction,  and  percolation  are  not  in¬ 
cluded  in  these  curves  because  they  are  site  specific.  The  effects  of  bottom 
friction  and  percolation  will  always  be  to  reduce  the  estimated  wave  height. 
These  curves  should  be  applied  only  in  areas  of  nearly  parallel  bottom  contours. 
Consequently,  refraction  will  also  only  reduce  wave  height. 

This  report  presents  a  method  for  estimating  the  significant  wave  height, 

Ha,  given  the  fetch  length,  F,  the  overwater  windspeed,  Uw  (see  U.S.  Army, 
Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1981),  and  the  water 
depth,  d,  neglecting  any  additional  wave  growth  in  shallow  water  due  to  the 
wind.  The  method  differs  from  two  recently  reported  methods— Seelig  (1980), 
who  presents  a  method  for  predicting  shallow-water  wave  height  given  deepwater 
wave  height,  H0,  peak  period,  Tp,  and  bottom  slope,  m,  and  Vincent  (1981), 
who  presents  a  method  for  calculating  the  depth-limited  significant  wave  height 
based  on  knowledge  of  the  deepwater  wave  spectrum — but  it  does  not  supersede 
the  use  of  these  other  two  methods.  The  report  provides  an  alternate  approach 
to  the  problem  of  shallow-water  wave  estimation  given  the  four  quantities  men¬ 
tioned  above. 

II.  WAVE  HEIGHT  PREDICTION  CURVES 

A  series  of  JONSWAP  spectra  was  generated  numerically  in  deepwater  condi¬ 
tions  for  varying  windspeeds  and  fetch  length,  and  propagated  into  shallow 
water  over  parallel  bottom  contours.  A  frequency-by-frequency  calculation  was 
made  at  various  depths  shoreward  applying  independently  the  wave  steepness 
limitation  criterion  (Kitaigorodskii,  Krasitskii,  and  Zaslavskii,  1975)  and  a 
shoaling  coefficient  to  each  spectral  component.  If  the  shoaled  wave  energy 
exceeded  the  limiting  value,  the  limiting  value  was  retained.  A  detailed  ex¬ 
planation  of  the  methodology  Involved  in  this  computation  is  presented  in  the 
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Appendix.  Resulting  spectra  at  grau-ally  decreasing  depths  for  a  given  case 
are  shown  in  Figure  1.  This  analysis  provides  the  wave  height  prediction 
curves  shown  in  Figure  2.  These  curves  provide  the  nearshore  significant  wave 
height,  Hs,  at  a  given  water  depth  which  is  related  to  the  total  energy,  E^, 
in  the  nearshore  wave  spectrum  by 


Ha  -  4Vet 

given  the  fetch  length,  the  overwater  windspeed,  and  the  deepwater  wave  height. 
Note  that  in  Figure  1  there  is  a  slight  shift  in  the  wave  period  toward  lower 
frequencies  as  the  spectrum  moves  into  shallow  water.  Later  work  will  attempt 
to  quantify  this  shift  and  incorporate  bottom  friction  effects. 

III.  USE  OF  CURVES 

n 

There  are  certain  restraints  on  the  use  of  the  curves  which  are  &b  follows: 

(1)  Curves  are  designed  to  be  used  for  fetch-limited,  wind-generated 
waves  in  deep  water  over  short  fetches,  i.e.,  up  to  62  miles  (100  kilo¬ 
meters)  . 

(2)  This  analysis  includes  only  the  wave  steepness  criterion  and 
shoaling.  It  does  not  reflect  other  energy  losses  such  as  refraction, 
friction,  or  percolation  (parallel  bottom  contours  are  assumed). 

(3)  The  fetch  length,  F,  is  strictly  the  straight-line  fetch 
unless  the  water  body  is  irregularly  shaped  where  the  fetch  would  be 
based  on  an  average  over  a  24°  quadrant. 


Figure  1.  Transformation  of  JONSWAP  spectrum  in  shallow  water. 
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Figure  2.  Dimensionless  fetch  versus  dimensionless 
wave  height  as  a  function  of  d/H0. 

(A)  To  calculate  the  adjusted  windspeed,  U^,  the  following  pro¬ 
cedure  should  be  used: 

(a)  If  windspeed  is  observed  at  any  level  other  than  33  feet 
(10  meters)  windspeed  on  land  or  water,  the  adjustment  to  the  33- 
foot  level  is  approximated  by: 


where  Uiq  is  the  10-meter  windspeed  in  meters  per  second,  Z 
the  height  of  wind  measurement  above  the  surface  in  meters,  and 
Uz  the  measured  windspeed  in  meters  per  second.  This  method 
is  valid  up  to  about  Z  ■  66  feet  (20  meters) .  If  the  windspeed 
was  measured  at  33  feet,  Uio  ■  U2. 

(b)  If  windspeed  was  measured  overland,  correct  to  overwater 
windspeed  by 
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Uw  -  1  •lUio  for  F  <  10  miles  (16  kilometers) 
Uw  ■  RUio  for  F  >  10  miles 


where  Uw  Is  the  overwater  wind speed  In  meters  per  second;  R 
Is  given  In  Figure  3.  If  wlndspeed  was  measured  overwater  and 
adjusted  to  a  10-meter  height,  ■  Uiq. 


Figure  3.  Ratio,  R,  of  wlndspeed  overwater,  1%,  to  wlndspeed 
overland,  Ul,  as  a  function  of  wlndspeed  overland, 

Ul  (after  Resio  and  Vincent,  1976). 


(c)  To  correct  for  wind  instabilities  over  fetch  lengths 
greater  than  10  miles; 


Ua  -  0.71  Ul-*3 

where  is  the  adjusted  wlndspeed  in  meters  per  second.  If 

the  F  s  10  miles,  U k  “  uw. 

(d)  To  correct  for  air-sea  temperature  differences, 


Ua  ■  Rt  for  F  >  10  miles 
UA  -  Ua  for  F  i  10  miles 

where  Ua  is  the  new  wlndspeed  adjusted  for  the  temperature  dif¬ 
ference;  Rij  is  given  in  Figure  4. 
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At  a  depth  of  3  meters ,  d/Ho  *  1.47,  providing  an  H  *  0.025  or  Hs  «  1.02 
meters.  The  peak  period,  Tp,  and  the  local  wavelength  would  increase  over 
that  at  a  7-meter  depth  and  currently  must  be  calculated  by  the  tables  given 
in  Appendix  C  of  the  Shore  Protection  Manual  (U.S.  Army,  Corps  of  Engineers, 
Coastal  Engineering  Research  Center,  1977). 

***************  LXAMPLE  PROBLEM  2 *************** 


GIVE!?;  The  wind  direction  is  predominantly  from  the  southwest  over  the  deep, 
irregularly  shaped  water  body  shown  in  Figure  5 .  The  windspeed  to  be  con¬ 
sidered  is  49.2  feet  (15  meters)  per  second  measured  on  top  of  an  instru¬ 
ment  shack  at  13  feet  (4  meters)  from  the  ground.  The  air  temperature  when 
these  conditions  occur  is  50°  Fahrenheit  (10°  Celsius)  and  the  water  tempera' 
ture  is  60°  Fahrenheit  (16°  Celsius). 


Figure  5.  The  fetch  length  for  this  irregularly  shaped  water  body  in  the  wind 
direction  is  determined  by  drawing  nine  radials  at  3°  increments 
centered  on  the  wir.d  direction  and  arithmetically  averaging  the 
radial  lengths  as  illustrated.  The  average  fetch  in  this  example 
is  approximately  22.2  miles  (36  kilometers). 


12 


FIND ;  The  significant  wave  height  at  a  16.4-foot  (5  meter)  depth  just  off  the 
coast  near  the  anemometer  site. 


SOLUTION :  The  fetch  is  found  by  averaging  over  a  24°  quadrant  since  the  body 
of  water  is  irregularly  shaped.  As  shown  in  Figure  5,  nine  radials  are 
constructed  at  3°  increments  and  the  average  fetch  length  of  22  miles  (36 
kilometers)  is  found. 

The  adjusted  windspced  is  found  following  the  steps  outlined  previously: 


(a)  Adjust  wind  from  the  4-meter  to  the  10-meter  level 
/10X1/7  /iox1/7 

Uio  *  (  2~)  UZ  "  (4^)  (15)  ■  17.1  meters  per  second 


(b)  Adjust  overland  wind  to  overwater  wind  with  R  from  Figure  3 

Uw  ■  RU10  ■  1.25(17.1)  -  21.4  meters  per  second 

(c)  Correct  wind  for  instabilities 

■  0.71  Uw1-23  *  0.71(21.4) 1 ‘ 23  ■  30.7  meters  per  second 

(d)  Correct  for  air-sea  temperature  difference  with  Rf  from  Figure  4 

Ua  ■  Rt  Ua  ■  1.17(30.7)  -  35.9  meters  per  second 

The  dimensionless  fetch,  x,  is 

-  _  8F  _  (9.8  m/s2) (36.000  m)  _  2?3  ? 
u£  (35.9  m/s)2 


The  deepwater  significant  wave  height  and  peak  period  are 

Hq  -  1.6  x  10- 


1-3  [36,000 
\  9.8 


g—  (35.9  m/s)  ■  3.5  meters 


T  .  35.9(273.7) 
P  3. 5(9. 8) 

At  a  5-  t'eter  depth 


J/3 


6.80  seconds 


d  5  .  ... 

H0  "  3.5  “  1,43 

In  Figure  2  at  x  •  273.7  and  d/Ho  -  1.43 
-  gHs 

H  -  -  0.012 

»* 


and 


Hb  - 


“  UA  (0.12) (35. 9)2 


9.8 


1.58  meters 
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APPENDIX 


METHODOLOGY  AND  GOVERNING  SPECTRAL  EQUATIONS 
1.  Deepwater  Representation  of  Fetch-Limited  Wav  Spectrum. 

A  spectrum  of  wind  waves,  generated  in  deep  water  for  a  long  period  of 
time,  is  limited  by  the  length  of  the  fetch  over  which  the  wind  Is  blowing. 
The  wind  will  generate  a  spectrum  with  a  shape  which  has  been  parameterlsed 
by  Hasselmann,  et  al.  (1973).  The  parameterization,  or  JONSWAP  spectrum,  pro 
vldes  a  functional  relationship  between  energy  and  frequency  as  well  as  the 
wind speed,  fetch  length,  and  width  of  the  spectral  peak: 


and 


E(f) 


f“5 


exp 


(f-fm)2 


(A-l) 


a 


era  for  £  <  £m 
ob  for  f  >  fm 


where 

E  *  energy  density 

F  ■  fetch  length 

f  *  frequency  of  wave  component 

fjn  *  frequency  of  spectral  peak  -  3.5g/(uio  x 

g  ■  acceleration  due  -to  gravity 

Ua  *  adjusted  10-meter  windspeed 

x  ■  dimensionless  fetch  -  gF/U^ 

a  ■  Phillips  equilibrium  constant  ■  0,076  x  “0.22 

Y  ■  ratio  of  maximal  spectral  energy  to  the  maximum  of  the  corresponding 

Pierson-Moskowitz  spectrum  *3.3 

oa  -  left-side  width  of  the  spectral  peak  -  0.07 

ob  -  right-side  width  of  the  spectral  peak  *0.09 


This  equation  provides  a  wave  spectrum  as  shown  in  curve  1  (Fig.  1),  with  a 
total  energy  equal  to  the  deepwater  significant  wave  height,  squared  over  16. 
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Reduction  In  Shallow  Wate\ . 


As  an  irvegular  wave  train  enters  transitional  and  shallow-vater  depths, 
the  presence  of  the  sea  bottom  causes  changes  In  wave  steepness  which,  due  to 
the  limitation  cn  wave  steepness,  lead  to  a  loss  of  wave  energy.  Kitalgorodskll, 
Krasltskil,  and  Zaslavskli  (1975)  suggest  that  an  upper  limit  of  energy  exists 
at  a  given  frequency  which  Is  a  function  of  depth  and  a: 


E(f)  -  ag2  f-5  (2ir)-4 

* 


(A-2) 


where 

Ch  tanh  (u>h  C^)  *  1 

d  ■  water  depth 

a  -  0.0081 

<t>  -  Ch{l  +  [2uh  Ch/sinh(2wh  Ch>  ] } 

a>k  “  2uf  /d/g 


This  equation  represents  a  stability  limit  or  "limiting  form  criterion"  on  a 
wave  component.  Kitalgorodskll,  Krasltskil,  and  Zaslavskli  used  a  value  of 
n  of  0.0081  based  on  field  data.  Recent  work  at  the  U.S,  Army  Engineer 
Waterways  Experiment  Station  (WF.S)  has  indicated  that  another  mechanism,  non¬ 
linear  wave-wave  interaction,  has  an  equivalent  effect  but  that  a  would  vary 
with  dimensionless  fetch  (gF/Uj[) .  The  application  of  this  theory  is  further 
v  *tlined  by  Vincent  (1981). 


Shoaling  of  a  wave  in  shallow  water  also  changes  wave  energy.  A  shoaling 
coefficient  can  be  calculated  as  in  the  Shore  Protection  Manual  (App.  C  in 
U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center,  1977)  for 
each  frequency  component  according  to  linear  theory: 


MO 


4TTd/L(f) 
sinh  4ird/L(f). 


(A-3) 


and  can  be  multiplied  by  the  deepwater  energy  at  each  frequency  band  to  obtain 
a  "shoaled"  spectrum, 


E(f)  shoaled  -  Ks(f)  E(f)  deep 


(A-4) 


3.  Determination  of  Shallow-Water  Energy  Spectrum. 


Figure  A-l  is  a  flow  chart  describing  the  solution  process  used  in  produc¬ 
ing  the  design  curves  presented  in  this  paper. 
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Ganarata  daapvattr  JOMSUAP  apac- 
truw  uaing  aaiiatlon  (A-l)  baaad 
on  vindapaad  and  fatch  langth. 


Figure  A- 


.  Flow  chart  Illustrating  the  use  of  equations  (A-l)  to 
(A-4)  in  generating  the  curves  presented  in  Figure  2. 
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